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We devise a theoretical description for the response of nitrogen molecules (N2) to ultrashort and intense 
x rays from the free electron laser (FEL) Linac Coherent Light Source (LCLS). We set out from a rate- 
equation description for the x-ray absorption by a nitrogen atom. The equations are formulated using all 
one- x-ray-photon absorption cross sections and the Auger and radiative decay widths of multiply-ionized 
nitrogen atoms. Cross sections are obtained with a one-electron theory and decay widths are determined 
from ab initio computations using the Dirac-Hartree-Slater (DHS) method. We also calculate all binding and 
transition energies of nitrogen atoms in all charge states with the DHS method as the difference of two self- 
consistent field calculations (ASCF method). To describe the interaction with N2, a detailed investigation 
of intense x-ray-induced ionization and molecular fragmentation are carried out. As a figure of merit, we 
calculate ion yields and the average charge state measured in recent experiments at the LCLS. We use a series 
of phenomenological models of increasing sophistication to unravel the mechanisms of the interaction of x rays 
with N2: a single atom, a symmetric-sharing model, and a fragmentation- matrix model are developed. The 
role of the formation and decay of single and double core holes, the metastable states of N 2+ , and molecular 
fragmentation are explained. 

PACS numbers: 33.80.-b, 33.80.Eh, 32.80.Aa, 41.60. Cr 
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I. INTRODUCTION 

The theory of the interaction of ultrashort and intense 
x rays with atoms and molecules has recently come into 
focus by the newly-built x-ray free electron laser (FEL) 
Linac Coherent Light Source (LCLS) in Menlo Park, Cal- 
ifornia, USAfii 2 - The first experiments clearly revealed 
how important a detailed understanding of the interac- 
tion of x rays with single atoms and small molecules is 
for their interpretation!^— To obtain a clear theoretical 
view on the physics of the interaction with such compar- 
atively simple systems represents a crucial piece of basic 
research; it is essential to form a foundation for studying 



the interaction of more complex systems like condensed 
matter or biomolecules with intense x rays. Due to the 
novelty of x-ray FELs, there is currently little informa- 
tion available. Naively, one may suspect that the impact 
of x rays on matter depends only on the absorbed dose, 
i.e., the energy deposited per unit mass of medium. This 
would mean that the detailed pulse characteristics, i.e., 
duration, spatial profile, and temporal shape, are irrele- 
vant. However, first theoretical and experimental studies 
of neon atoms give a lucid explanation why this is too 
simple a picturej 3 -^^ In this work, we report on the the- 
oretical foundation of the interaction of intense x rays 
with small molecules and particularly emphasize the ad- 
ditional challenges one faces over describing interactions 
with single atoms due to molecular fragmentation. 



a ' Corresponding author; Electronic mail: christian. buth@ web. del 



One of the simplest molecules, nitrogen (N2), was stud- 
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FIG. 1. (Color online) Classification of key types of core holes 
in N2 in the localized core hole picture! 12 ' 13 (a) a singly- 
occupied spatial core orbital on one of the atoms of N2: a 
single core hole (SCH), (b) a fully-depleted spatial core or- 
bital on one N atom: a double core hole (DCH) on a single 
site (ssDCH), and (c) a singly-occupied spatial core orbital 
on both N atoms: a DCH on two sites (tsDCH). 



ied extensively in the first experiments at LCLS^— ^ and 
the multiphoton nature of the interaction was charac- 
terized. Specifically, Hoener et alA examined the in- 
teraction of x-ray pulses of constant energy but with 
varying duration and found in the ion yields the effect 
of frustrated x-ray absorption. Next, Fang et al& in- 
vestigated the formation of (multiple) core holes; the 
most important core-hole types were found to be sin- 
gle core holes (SCHs), double core holes (DCHs) on a 
single site (ssDCH) and on two sites (tsDCH) ; 12 i 13 all of 
them are depicted in Fig. [TJ The tsDCH is a novel type 
of DCH which can be produced at third-generation syn- 
chrotrons only with very low probability its abundant 
production has only become feasible with the emerging 
x-ray FELs^&iS^ Further, Cryan et al£^ investigated 
the formation and decay of DCHs in laser aligned N2 and 
the energy shift of DCHs with respect to the main SCH 
line where the synchronization between the optical laser 
and LCLS was achieved with the method of Glownia et 
al.^ they found a clear sign of a nonspherical angular 
distribution of the Auger electrons in the molecular ref- 
erence frame. Specifically, the experiments of Fang et al. 
and Cryan et al. open up novel perspectives for chemical 
analysis i 12 ' 13 i 15 ' 16 Finally, Hau-Riege et al& investigated 
near- ultraviolet luminescence of N2 . 

The prototypical N2 molecule has also been studied ex- 
tensively with other x-ray/xuv sources. Low- flux x rays 
from third-generation synchrotrons have been used to 
study the inner-shell properties of N ^ 17 ' 18 and the dis- 
sociative one-x-ray-photon absorption cross section was 
measured^ 7 - The production of SCHs by x-ray absorp- 
tion results in Auger electron spectra from SCH decay 
that were classified in Rcfs. [l9l - l23l . For photon ener- 
gies below the SCH-formation threshold, we note re- 
cent investigations at the Free Electron Laser in Ham- 
burg (FLASH) The study of intense xuv absorption 
by N2 at FLASH has a different character from x-ray ab- 
sorption by N2 at LCLS as the former radiation primarily 
targets valence electrons whereas the latter almost exclu- 
sively interacts with core electrons. 

In contrast to the one-photon physics at present syn- 
chrotrons, the interaction of the x rays from FELs with 



molecules is nonlinear. Namely, for intense x rays, sev- 
eral one-photon absorption events are linked by the de- 
cay widths of the created holes. So not only two or more 
photons are absorbed but the result of the absorption 
also depends on the time-difference between the individ- 
ual one-photon absorption events. If it did not depend 
on the time-difference, the outcome of multiphoton ab- 
sorption would only depend on the number of absorbed 
photons and not on the timing of the processes, i.e., the 
processes needed not be considered as an entity; how- 
ever, due to the finite lifetime of inner-shell holes, the 
two events are linked and cannot be treated separately. 

In this paper, we analyze theoretically the interac- 
tion of intense x rays with N2. They represent the next 
step from studying the impact of x rays on single atoms 
like neonj££iii Molecules in intense x-ray FEL light of 
high enough photon energy can be multiply core-ionized 
by sequentially absorbing several x-ray photons. Direct 
valence ionization occurs only rarely because the asso- 
ciated cross sections are orders of magnitude smaller 
than cross sections of core electrons <2£r— In very in- 
tense x rays, this may even lead to a hollowing out of 
molecules, i.e., core shells may be completely depleted 
due to ssDCH formation. For depleted shells, further ion- 
ization is drastically reduced and can even lead to trans- 
parency for x raysi^ 9 . Significant ionization of hollowed 
out molecules resumes only when Auger decay refills the 
core shells. We will theoretically analyze the copious- 
ness and influence of core holes for N2 in intense x rays. 
We set out with an investigation of isolated N atoms 
and compute ion yields and the average charge state. 
Then, we consider N2 molecules for which we have an- 
other possibility that has significant impact on observ- 
ables; namely, as long as the molecular ion can be re- 
garded as an entity, i.e., dissociation has not yet sepa- 
rated the molecule into fragments, the valence electrons 
are shared by both N atoms. This leads to a substantial 
redistribution of charge in decay processes of core-ionized 
molecule s 17 ' 19 " — which has a significant influence on the 
observed ion yields from N2 (and many other) quanti- 
ties. Dicationic states N 2+ are predominantly populated 
by Auger decay of SCHs; they are frequently metastablc 
with respect to dissociation with comparatively long lifc- 
times i 23 i 30 In this case, dissociation on a rapid timescale 
only resumes when more electrons are removed from the 
dication. This interdependence of multiple x-ray absorp- 
tion and nuclear dynamics — which eventually leads to 
molecular fragmentation — poses a challenge to theoret- 
ical molecular physics. In the present study, we use a 
series of phcnomenological models to unravel the compli- 
cated electronic and nuclear processes; it complements 
the phcnomenological molecular rate-equation model 31 
described in Ref. 

The paper is structured as follows: we present a the- 
oretical study of the absorption by N2 of ultrashort and 
intense x-ray pulses from LCLS at a wavelength of 1.1 nm 
(1100 eV photon energy) with a nominal pulse energy 
of 0.15 mJ for FWHM pulse durations of 4 fs and a pulse 
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energy of 0.26 mJ for FWHM pulse durations of 7, 80, 
and 280 fs. For multiply-ionized N atoms, we use ab ini- 
tio computations to determine Auger and fluorescence 
decay widths and the corresponding transition energies, 
electron binding energies, and one- x-ray-photon absorp- 
tion cross sections [Sec. HI]. We formulate a rate-equation 
model for an N atom and study the sequential absorp- 
tion of multiple x-ray photons in Sec. IIII1 Ion yields of 
a single N atom are predicted and show a strong devia- 
tion from the experimentally measured ion yields of N2. 
The atomic rate-equations are, therefore, extended phc- 
nomenologically in Sec. IIVI to treat molecular fragmen- 
tation which is the cause for the observed differences be- 
tween the ion yields of atoms and of molecules. For this 
purpose, a series of models of increasing sophistication 
is developed. In Sec. [Vj we explain the phenomenology 
of N2 in intense x rays in terms of elemental molecular 
processes. Conclusions are drawn in Sec. IVII 

Our equations are formulated in atomic units^ For 
the conversion of a decay width T in electronvolts to a 
lifetime r in femtoseconds, we use the relation r = |l = 
0.658212 cVfe ^ Au dctails of the m odels and calculations of 
this paper are provided in the supplementary materials^ 



II. ENERGIES AND DECAY WIDTHS OF 
MULTIPLY-IONIZED NITROGEN ATOMS 

In this section, we study the energy levels of multiply- 
ionized N atoms and their decay widths and transition 
energies of Auger and radiative decay channels. They 
will be used to describe the interaction of atoms with 
x rays in rate-equation approximation in Sec. IIIII From 
the energy levels of multiply-charged N atoms, we iden- 
tify the energetically accessible ionization channels for 
a specific x-ray photon energy because only those elec- 
trons are available for direct one- x-ray-photon ionization 
that have a lower binding energy than the photon energy. 
The transition energies for Auger decay and x-ray fluo- 
rescence allow an identification of the relevant processes 
in experimental electron and photon spectra. For ele- 
ments from the second row of the periodic table, like 
nitrogen, only Auger decay is important, if a competing 
Auger decay channel exists in addition to x-ray fluores- 
cence [Tables IV, VI and IX, XI in Ref. [H|. 

Energy levels, transition energies, and decay widths for 
multiply-ionized N atoms are calculated using the Dirac- 
Hartree-Slater (DHS) method i 33 i 34 In the energy calcu- 
lations, we use the DHS wave functions from the appro- 
priate defect configurations to compute the expectation 
values of the total Hamiltonian including the Breit inter- 
action^ i.e., a first-order correction to the local approx- 
imation is made. We also include the effects of relaxation 
by performing separate self-consistent field (SCF) calcu- 
lations for the initial and final states [ASCF method] and 
take the difference of the total energies to obtain bind- 
ing energie s 33 ' 35 or transition energies, respectively. We 
calculate Auger and radiative decay widths with pertur- 



bation theor y 3 i 36 from frozen-orbital and active-electron 
approximations as opposed to the relaxed-orbital calcu- 
lations used to determine energies. The Auger matrix 
elements are calculated with DHS wave functions for the 
initial state with a proper defect configuration. Contin- 
uum wave functions are obtained by solving the Dirac- 
Slater equations with the same atomic potential as for the 
initial state. With this treatment, the orthogonality of 
the wave functions is assured. The radiative decay widths 
are calculated in multipole expansion using the relativis- 
tic length gauge for the electric multipoles^ The radial 
matrix elements are also calculated using DHS wave func- 
tions from the initial state. The effects of relaxation on 
the total ET-shell Auger and radiative decay widths are 
less than 10 % for atoms with a AT-shell vacancy^ For 
a few times ionized N atoms, the effects are even smaller 
of the order of only a few percent. 

The rclativistic calculations of energies and decay 
widths for a N atom yield values in the jj coupling 
scheme, i.e., they are specified with respect to electronic 
configurations which are labeled by ls e 2s m 2p 1 ^ 2 ■ 
Here, < £,m,n 1 / 2 < 2, and < n 3 / 2 < 3 are 
the occupation numbers of the respective orbitals where 
the subscripts 1/2 and 3/2 denote the value of the to- 
tal angular momentum quantum number j£^ By bas- 
ing our method on DHS wave functions, we develop a 
more general approach than previous investigations, e.g., 
Ref. [33, by treating the atomic electronic structure in 
one-electron approximation including rclativistic effects. 
The DHS approach yields good accuracy throughout the 
periodic table and facilitates to assess, and if necessary 
incorporate, such effects. For light elements like nitrogen, 
rclativistic effects are small and have only minor impor- 
tance for our modeling. Specifically, spin-orbit coupling 
is small which allows us to determine energies and decay 
widths in LS coupling — where the electronic configura- 
tions are denoted by Is 1 2s m 2p n with n = rix/2 + ^3/2 — 
from the corresponding quantities in jj coupling scheme. 
Scalar rclativistic effects arc, nonetheless, accounted for 
in this procedure. 

In the supplementary materials, 28 we discuss dctails of 
the probabilistic averaging of energies and decay widths 
in the jj coupling scheme and we provide following data 
for a multiply-ionized nitrogen atom: the electron bind- 
ing energies are gathered in Tables III, VIII, and XIII; we 
list the Auger decay widths in Tables IV and IX with the 
corresponding transition energies provided in Tables V 
and X; the x-ray fluorescence decay widths are presented 
in Tables I, VI, and XI with the respective transition en- 
ergies collected in Tables II, VII, and XII. Our Auger 
decay widths for SCHs in Table IV are in reasonable 
agreement with the early computations^ of Bhalla4£ 



III. ION YIELDS OF A NITROGEN ATOM 

To describe the time-evolution of the absorption of 
x rays and the induced decay processes by a single 
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N atom, we use rate equations. This approach has proven 
useful in many problems of light-matter interaction and 
also for the interaction of x rays from free electron lasers 
with atomSi 11 ! 41 " — The justifications for this approxima- 
tion are, first, that the envelope of the x-ray pulse varies 
slowly with respect to the change of populations of states, 
i.e., ionization and decay processes happen almost in- 
stantaneously on this time scale and, second, coherence 
effects are small; the x rays induce essentially only one- 
photon absorption processes which are, however, linked 
by inner-shell hole lifetimes. 

The rate equations for a single nitrogen atom 2 ^ com- 
prise only electron removal terms and no electron addi- 
tion terms because only the former are relevant for our 
study and the considered nitrogen charge states are al- 
ways neutral or cationic. We include all possible one- 
x-ray-photon absorption processes in the independent 
electron approximation; the cross sections are obtained 
in one-electron approximation with the Los Alamos Na- 
tional Laboratory Atomic Physics Codes ^^ 45 All Auger 
decay widths of Sec. [TT] are used, however, x-ray fluo- 
rescence is only accounted for when there is no compet- 
ing Auger channel because Auger widths are typically 
orders of magnitude larger than the corresponding flu- 
orescence widths^ Other contributions are neglected, 
specifically, cross sections for the simultaneous absorp- 
tion of two x rays are tiny^ Likewise more complicated 
and less important processes like shake ofl& and double 
Auger decay^ are not incorporated. 

The probability to find an N atom at time t in an elec- 
tronic configuration ls l 2s m 2p n is denoted by Pimnif) 
where i — m = 2, n = 3 stands for the neutral atom. 
The x-ray photon flux at time t is given by Jx(t)i i.e., 
the number of x-ray photons that pass a unit area per- 
pendicular to the beam in unit time^- The cross sec- 
tion at a chosen x-ray photon energy for the ionization 
of an electron of an N atom with electronic configura- 
tion Is 1 2s m 2p n is denoted as ai mn with < £, m < 
2 A < n < 3. With these quantities, we formulate 
the rate equation to describe the depletion of the ground 
state of nitrogen due to x-ray absorption: 

=-<T22sP22 3 (t)J X (t). (1) 

The total (Auger plus x-ray fluorescence) decay width of 
a singly core-ionized state Is 1 2s 2 2p 3 is specified as Fi23. 
The Is ionization cross section of an N atom is denoted 
by CTi23<_223- Using these parameters, we find for the 
single core-hole rate 



dP 123 (t) 
dt 



= (o"l23^223 P 2 23{t) ~ °~\23 Pl23(t)) Jx{t) 



Tl23 Pl23(t) 



and analogously for the double core-hole rate 
dPo 23 (i) 



(2) 



dt 



C023<-123 Pl23{t) — C023 -P023W) Jx(t) 
- To23 Po23(t) . 



The three rate equations (J]), ©, and (j3)) already are 
sufficient to describe the formation of SCHs and DCHs 
in an N atom. They form a system of ordinary linear 
differential equations which is solved numerically with 
Mathematical for the following initial conditions: before 
the x-ray pulse, the N atom is in its electronic ground 
state, i.e., P 2 23(~oo) = 1, and all other charge states 
are unpopulated, i.e., P[ mn {—oo) = for ->(£ = m = 
2 A n = 3). Many more rate equations — 36 in total — 
are needed to describe the ionization and various decay 
channels for all charge states from a neutral to a septuply- 
ionized N atom. 

With the rate-equation description at hand, we can 
compute obscrvables. The probability to find an initially 
neutral N atom in a cationic state with charge < j < 7 
after the x-ray pulse is over (t — > oo) is 



Pj = 



E 



Plmn (oo) 



(4) 



£+m+n—7—j 
0<f,ra<2AO<n<3 



We define atomic ion yields to be normalized charge-state 
probabilities with 1 < j < 7 by 



Y; 



Pi 



1-Pn 



(5) 



where 1 — Pq = ^2 Pj is the probability to create a cation. 

k=l 

They are experimentally accessible and are a meeting 
point with theoryj^ The Yj can be used to determine an 
average charge state q via 
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(0) 



(3) 



It is a useful gross quantity to characterize the x-ray in- 
duced ionization and quantifies the average amount of 
radiation damage caused^ 

Exemplary charge-state probabilities (U]) for an iso- 
lated N atom irradiated with intense x rays are shown in 
Fig. [5^ for a range of photon energies from to 2000 eV£& 
At low photon energies, certain ionization channels and, 
therefore, higher charge states are not accessible because 
the photon energy is not large enough to ionize inner 
electronic shells. Specifically, energies of shells move to 
higher binding energies with increasing charge state of 
the ions such that they move out of reach for a one- 
photon ionization process with fixed x-ray photon en- 
ergyAi This is particularly clearly discernible in the range 
from 400 to 700 eV where steps indicate openings of ion- 
ization channels. Above 1200 cV, higher charge states 
decrease monotonically and the ground-state remains 
more and more populated because the cross sections have 
dropped substantially in magnitude such that the x rays 
do not ionize as effectively as before anymore. 

The average charge state q [Eq. ©] for a single N atom 
in intense x rays is plotted in Fig.[5}3 for the same range of 
photon energies as in Fig. [5^,. The above behavior of the 
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FIG. 2. (Color) An N atom after irradiation by x rays of constant spatial beam profile with a photon fluence (photon flux 
integrated over all times) of 8 x 10 11 /mi -2 and a Gaussian temporal shape (O with FWHM duration of 2fs at varying 
photon energy uix- (a) Probabilities Pj(tox) from Eq. Q to find an N atom in specific charge states: neutral (black), singly 
(red), doubly (green), triply (blue), quadruply (magenta), quintuply (cyan), sextuply (purple), and septuply (orange), 
(b) Average charge state q(u)x) from Eq. ©. 



charge-state probabilities Q is mirrored in the depen- 
dence of q on the x-ray photon energy. For low photon 
energies, q rises quickly due to large cross sections but its 
maximum value is limited by the energetic accessibility of 
electronic shells; it peaks, first, around 100 eV when the 
valence shells of N atoms can be fully depleted (Table III 
in Ref. [28h producing quintuply ionized final states and 
thus q — 5. For higher photon energies q drops again due 
to decreasing cross sections. Starting at about 400 eV 
also core electrons of nitrogen can be accessed leading to 
a second peak at about 700 eV when all electrons can be 
ionized. This peak does not reach the maximum value 
of q = 7 because the x-ray fluence is too low to fully strip 
an N atom of its electrons. This is in contrast to the 
first peak where the maximum value is actually reached. 
Above 700 cV photon energy, q drops monotonically due 
to decreasing cross sections. 

To describe the interaction of finite pulses of limited 
transversal extend — along the x and y coordinates with 
respect to the propagation direction along the z axis — 
with an ensemble of atoms, we need to integrate over the 
temporal and spatial profiles of the incident x-ray pulse 
in the focal region. The LCLS is a FEL which operates 
by the self-amplification of spontaneous emission (SASE) 
principle^— which produces chaotic radiation pulses 
with random spikes. In principle, the pulse structure of 
SASE pulses needs to be taken into account in nonlinear 
optical processes. However, as exemplified for atoms in 
Ref. [ll|, the impact of chaotic light on linked one-photon 
processes — which are the predominant mode of interac- 
tion of intense x rays with nitrogen — has little impact. 
Therefore, there is no need to treat the explicit form of 
SASE pulses in this problem [see also Table [TT| for a com- 
parison of results for a Gaussian pulse with results from 
SASE pulses]. Instead, it is sufficient to use a Gaussian 
temporal profile for the rate of photons 

r x W = r x ,oc- 41n2( ^ )2 , (7) 

with a peak rate of Tx,o = 2 J ^ ^ for n p h photons 



in the pulse and a full width at half maximum (FWHM) 
duration of rx- The number of photons ?i p h is linked to 
the total energy in the pulse Ep via n p h = for approx- 
imately monochromatic x rays with photon energy wx- 
The LCLS beam is modeled as a Gaussian bearn^ with a 
long Raylcigh range which implies that the spatial profile 
of the beam along its axis, the z axis, does not vary appre- 
ciably over the acceptance volume of the ion spectrome- 
ters. Therefore, the x-ray flux is approximately constant 
along the z axis and thus we may disregard its z depen- 
dence in what follows. There is, however, a substantial 
variation of the x-ray flux transversally to the beam axis. 
The transversal spatial profile of LCLS pulses was de- 
rived from ablation of solids in the FEL beam^i - — that 
yielded an elliptic profile with a FWHM length of the 
major and minor axes of g» x = 2.2 /mi by g y = 1.2 /mi, 
respectively. We have the following analytic expression 
for the transverse beam profile 

it is normalized to unity with respect to an integration 

oo oo 

over the entire x—y plane, i.e., J J g(x, y) Ax dy = 1. 

— oo — oo 

In our computations we find that an integration in the x— 
y plane over an area of 10 2 /tm 2 centered on the beam axis 
yields converged ratios of DCHs versus SCHs. Finally, 
the position-dependent x-ray flux follows from Eqs. ([7]) 
and ([8]) to 

J^(x,y,t) = g(x,y)T x (t) . (9) 

We use only a fraction of the nominal pulse energy^ to 
determine n p h to account for losses in the experimental 
setup. This percentage is determined in Sec. IIVBI by a 
fit of experimental data with our theoretical predictions. 
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(b) 



N + N 2+ N 3+ N 4+ N 5+ N 6+ N 7+ 



FIG. 3. (Color online) Ion yields of (a) an N atom and (b) an N2 molecule irradiated by x rays of a FWHM duration of 280 fs 
and a photon energy of 1100 eV. In (a) we show theoretical results from the single atom model of Sec. IHII for the x-ray flux ((9|. 
We assume that only 31% of the nominal pulse energy of 0.26 mj are actually available in the LCLS experiment [Table H]. In 
(b) we show experimental data for an N2 molecule. The comparison of atomic Y, with molecular Yj ion yields for nitrogen 
clearly demonstrates the impact of molecular processes during and after the interaction with the x rays. 



IV. ION YIELDS OF A NITROGEN 
MOLECULE 

In the previous Sec. IIII| we discussed the ion yields of 
an isolated N atom which we denoted by Yj for the posi- 
tive charge j £ {1, . . . , 7}. In this section ion yields from 
an N2 molecule — after fragmentation into atomic ions — 
arc calculated to which we reference by Yj . Figure [3] 
shows the theoretical ion yields for an N atom alongside 
the experimental data for an N2 molecule; we find a stark 
difference between both bar charts. The trend found in 
the ion yields for a specific set of parameters solidifies, 
if we vary the pulse duration [Fig. 0] where we use q 
[Eq. (J6])] as a figure of merit. 

At first sight, one may assume that the absorbed dose 3 - 
and thus the ionization of a sample of N2 is very simi- 
lar when the x-ray pulses have the same total energy 
and central wavelength and only one-photon absorption 
processes are included in the rate equations [Sec. HVj . 
If there was only a single one-photon process with cross 
section <7i, then the probability of ionization after the 

x-ray pulse would only depend on the x-ray fluence, i.e., 

00 

Ci / Jx(t) dt, and thus be obviously independent of the 

—00 

time evolution of the x-ray flux Jx(i)- In Fig. [4] we 
display q with respect to the x-ray pulse duration deter- 
mined from the experimental data and for an N atom. 
One may expect naively that the absorbed dose is simi- 
lar in all cases and q is the same for all pulse durations, 
i.e., a horizontal line in Fig. |4] However, the contrary is 
true; we find a substantial dependence of q on the pulse 
duration. This finding is indicative of (nonlinear) multi- 
x-ray-photon physics in the sample, i.e., the absorption 
of several x rays. 

As expected from Fig. f3J there is a large discrepancy 
between q for an N atom and the experimentally ob- 
served q for N2. These differences cannot be ascribed to 
molecular effects on photoionization cross sections^ be- 
cause they are small for the chosen photon energies which 



are high above the ionization thresholds of multiply- 
charged N 2 i 18 ' 55 ~ — Likewise fluorescence and Auger de- 
cay widths in molecules are usually well-represented by 
the respective values of the isolated atoms. 

The differences between the molecular ion yields and 
the atomic ion yields in Fig. [3] are caused by the inter- 
play of core ionization of N2 with Auger decay, and nu- 
clear dynamics. Namely, x-ray absorption by N2 predom- 
inantly leads to SCH formation and Nj" cations; unless 
more x rays are absorbed shortly after, the SCH Auger 
decays producing the molecular dication N 2 , -1- . The N?, + 
fragments as follows: 

N 2+ — > N + N 2+ (10a) 
N 2+ — > N+ + N + , (10b) 

where the probabilities for the channels (|10a[) and (|10b() 
are 0.26 and 0.74, respectively— The N 2 , 4 " are an extreme 
case in the sense that the impact of molecular fragmen- 
tation is very pronounced on the ratio of N + yield versus 
N 2+ yield. The existence of the two fragmentation chan- 
nels (|10|) offers a lucid explanation for the observed differ- 
ences in Fig. [3] between atomic and molecular ion yields. 
For low x-ray flux, essentially no further x rays are ab- 
sorbed prior to the Auger decay of a SCH in an N atom, 
i.e., the initially neutral atom ends up as N 2+ . How- 
ever, for Auger decay of a SCH in N2, the two charges, 
may be shared between both atoms (|10b[) . This fact ex- 
plains the large amount of N 2+ ion yield for the atom 
versus the large N + ion yield for the molecule in Fig. [3] 
A similar argument explains also the different ratio of 
the ion yields of N 3+ versus N 4+ ; namely, in the case 
that x-ray absorption is fairly slow — such that Auger de- 
cay and fragmentation mostly occur before another x ray 
is absorbed — a second x ray is more likely absorbed by 
a N + fragment than a N 2+ fragment as the former are 
more abundant than the latter. A subsequent Auger de- 
cay produces then preferably N 3+ instead of N 4+ . 

Despite the seemingly clear view on the physical rea- 
sons behind the deviations between Yj and Yj in Fig. [5J 
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one should bear in mind that posses metastable di- 
cationic states with respect to dissociation. The lifetime 
of N| + can extend even to a few microseconds^ for the 
dissociation through tunneling into the N + + N + con- 
tinuum. In a laser-dissociation experiment^ 80 fs and 
120 fs for the N+ + N+ and N + N 2+ fragmentation chan- 
nels, respectively, were measured which are on a compa- 
rable time scale with the LCLS pulse durations discussed 
here. Consequently, the impact of the metastable dica- 
tionic states is to slow down fragmentation with a time 
constant of ~100fs. This assessment is also supported 
by the fact that less than 2% of the N + peak in the time- 
of-flight spectra were ascribed to undissociated N| + in 
Ref. |j. For metastable states, dissociation on a rapid 
timescale may only resume when more electrons are re- 
moved from the dication. In the following Sees. IIV Al and 
IIV B[ we present several phcnomcnological models to ac- 
count for molecular fragmentation. 



A. Symmetric-sharing model 

Here we develop the simplest model of breakup of N2 
after irradiation by x rays in which charges produced by 
x-ray absorption on both N atoms are shared equally 
upon molecular fragmentation for an even number of 
charges; for an odd number of charges, one of the frag- 
ments carries one excess charge with respect to the other 
fragment. This model aims to contrast the extreme view 
of the single-atom model for ion yields of N2 from Sec. IIIII 
in which no redistribution of charge over the atoms in a 
molecule is accounted for in the course of the fragmen- 
tation process. Symmetric sharing is suggested in our 
situation by the dissociation of N^" 1 " [Eq. (|10p j. the anal- 
ysis of the potential energy curves of N| + which indi- 
cate a most probable dissociation into N + + N 2+ ,— and 
laser-dissociation experiments of Nj 4 " for which symmet- 
ric sharing was found to be the most prominent frag- 
mentation channel . 60 ' 61 As a redistribution of charge — 
that is produced treating the atoms independently — 
never increases the maximum charge found on the in- 
volved atoms, the single-atom ion yields of Sec. IIIII lead 
to an upper limit of the average charge state q from N2. 
Here we investigate the other extreme position which is 
to assume that the electrons of a charged N2 molecule 
are always shared equally between both N atoms upon 
fragmentation. In other words, the charge on the two 
atoms of N 2 is determined independently and then re- 
distributed between the two N atoms^ Such an equal 
sharing of charges leads to the smallest amount of pos- 
itive charge per atom and thus to the lowest q. Hence 
it represents a lower limit to the observed charge-state 
distribution. 

The q from the symmetric-sharing model is displayed 
in Fig. [4j Interestingly, we find that the lower limit is 
approached quite closely by the experimental data for 
short x ray pulses. However, the slow progression of the 
lower limit of q with the pulse duration leads to a signif- 



icant deviation for longer pulses. This indicates that, for 
short pulses, there is mostly equal sharing of the charges 
between the atoms in N2 while there is a different mech- 
anism at work for longer pulses. Of course, this is related 
to the fact that the importance of the various processes 
entering the atomic rate equations of Sec. lIIII are weighted 
differently depending on the actual LCLS pulse parame- 
ters. As will be discussed in detail in Sec. [V] the princi- 
pal reasons for the observed deviations is that during the 
long x-ray pulses, Auger decay processes mostly occur 
prior absorbing more x rays and fragmentation may hap- 
pen in between. However, for short pulses, the molecule 
charges up quickly and nuclear expansion and fragmen- 
tation are of little importance on the time scale of the 
shortest (4fs) pulse. In this case, only after the pulse 
is over, there is sufficient time for the induced nuclear 
wavepacket to evolve leading to molecular breakup. 



B. Fragmentation-matrix model 

To account for fragmentation of N2 in the course of the 
interaction with x rays in a more realistic way, we make 
the following heuristic ansatz. We use the probabilities 
from the atomic rate equations after the x-ray pulse is 
over to form probabilities to find N2 in various charge 
states treating the two N atoms as independent. The 
fragmentation of the molecular ion is then assumed to be 
determined only by the charges on the two atoms in terms 
of fixed ratios independent of the timing of the interac- 
tion with the x rays. This approximation is motivated 
by realizing that the fragmentation of N 2-1 " — for which 
we have the two channels of Eq. (|10|) — depends only on 
the amount of N| + produced in the course of the interac- 
tion with the x-ray pulse. A fixed fragmentation constant 
then accounts for the redistribution of charges. Intro- 
ducing constants for all relevant fragmentation patterns 
of N2 in all possible charge states leads to a fragmen- 
tation matrix which transforms the molecular probabili- 
ties from the independent atom treatment into molecular 
probabilities that account for charge redistribution from 
which ion yields arc determined that can be compared 
with experimental data. Using a fixed fragmentation pat- 
tern clearly implies that we assume that fragmentation 
always occurs in the same way for a given independent- 
atom molecular charge state after the x-ray pulse is over, 
i.e., we neglect all dependencies on the actually involved 
intermediate states which change for different x-ray pulse 
parameters and dissociation during the x-ray pulse and 
focus exclusively on the achieved final independent-atom 
molecular charge state. 

For a mathematical formulation of the fragmentation- 
matrix model, we set out from the atomic probabilities 
to find an N atom in various charge states after the 
x-ray pulse is over Pj with < j < 7. They are readily 
obtained from the solution of the atomic rate equations of 
Sec. IIIII via Eq. |QJ. Independent- atom molecular prob- 
abilities P[: of N2 for finding charge state i on the left 
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N atom and charge state j on the right N atom follow 
then immediately from P[- = Pi Pj = P^ for < i,j < 7; 

they are aggregated in the vector P' . The fragmentation 
matrix is denoted by J 1 ; it transforms the P' into molecu- 
lar probabilities that account for charge redistribution P 
via 

P = FP'. (11) 

The vector P is formed by the molecular probabili- 
ties P i:j = P r i to find N { 2 i+3)+ . Nota bene, Eq. §\A§ 
needs to be applied for a single x-ray intensity; volume 
averaging (|8|) is performed for the transformed result. 
The atomic probabilities Pj [Eq. (J4j) ] are recovered from 
Eq. (HH) by choosing F = 1. Summing over one index 
of Pij then yields 

7 

7 

because Pj = 1 holds. Relation (|12|) is a general prop- 

i=0 

erty of any P from a meaningful F to find the probabili- 
ties Pj of atomic fragments after molecular breakup. For 
an arbitrary F, we find by setting P' in Eq. (|lll) to Carte- 
sian unit vectors that the sum of the matrix elements in 
each column of F is unity, i.e., the sum of the elements of 
a vector remains the same (unity for probabilities) upon 
transformation with F. 

To find a form of F that introduces a redistribution 
of charge, we make the additional assumption that a 
molecular ion will dissociate into fragments with simi- 
lar charge states. Specifically, we assume that a molec- 
ular ion with an even number of electrons predomi- 
nantly fragments symmetrically and, for an odd num- 
ber of electrons, it fragments such that only one excess 
charge is found on one of the two fragments. As dis- 
cussed previously, the independent-atom view [Sec. IIII] 
and the symmetric-sharing view [Sec. IIV A| are two ex- 
treme positions with the former representing the up- 
per limit of the average charge state [Eq. ©] and the 
latter being its lower limit. Hence a redistribution of 
charge between the independent-atom channel and the 
symmetric-sharing channel, with a certain amount of 
probability shifted from the former to the latter, is able 
to represent the impact on q even if more fragmenta- 
tion channels are present for a specific molecular charge 
state. In other words, we realize that only the cou- 
pling of P[j for \i — j\ > 2 to more symmetric charges 
states needs to be considered for our purposes and a 
fraction < fs%j\ < 1 of the probability of the charge 
state {i,j} goes over into {k,l} = {(i + j)/2, (i + j')/2} 
for even i+j and into {k,l} = {(i+j -l)/2, (i+j + l)/2} 
for odd i + j.— This kind of redistribution of probabil- 
ity can be introduced easily using a transformation ma- 
trix T- '{i.j} which differs from the identity matrix by set- 
ting ( :F {ij})( tJ ),( 4 ,i) = {^{%3})(iA,(JA = 1 ~ f{ij} and 



r x [fs] 


E P [mJ] 


/{0,3} 


/{0,4} 


/{2,4} 


/{0,5} 


/{2,5} 


/{0,6} 


280 


0.31 x 0.26 





0.54 


1 








1 


80 


0.25 x 0.26 





0.61 


1 


0.26 


0.17 


1 


7 


0.16 x 0.26 





1 


1 


1 


1 


1 


4 


0.26 x 0.15 


1 


1 


0.48 


1 





0.65 



TABLE I. Fragmentation constants for the most probable 
channels for the four pulse durations. Here, /{o,2} = 0.74 in 
all cases [Eq. iflQ]!]. The LCLS FWHM pulse duration is r x 
[Eq. {7)], the LCLS photon energy is 1100 eV, and the actu- 
ally available pulse energy in the LCLS experiment is Ep. A 
nominal pulse energy of 0.15 mJ is specified for 4 fs pulses and 

0. 26 mJ for the remaining three pulse durations. 

(^{^■^(MMij) = (■ ?7 {i 1 j})a,fe),(j>-) = f{i,j}- Thc sct of 

21 molecular charge states {i, j} for which \i — j\ > 2 
holds is 

K = {{0, 2}, {0, 3}, {0, 4}, {0, 5}, {0, 6}, {0, 7}, 

{1,3},{1,4},{1,5},{1,6},{1,7},{2,4}, (13) 
{2, 5}, {2, 6}, {2, 7}, {3, 5}, {3, 6}, {3, 7}, 
{4, 6}, {4, 7}, {5, 7}}. 

For all k 6 K, we find a J~k from which we construct 
the fragmentation matrix via -F 1 = H Fk- As ZFk and 

fceif 

J-i commute for all k, I G K and matrix multiplication is 
associative, F does not depend on the sequence of mul- 
tiplication of the T^. The set K [Eq. (fTB"]) ] contains all 
necessary molecular charge states to redistribute charges 
as is done in the symmetric-sharing model of Sec. IIV AI 
Consequently, the probabilities of Sec. IIV Al are recovered 
by choosing /fc = 1 for all k S K and using Eq. (|12[) . 

Having asserted us that the fragmentation-matrix 
model correctly integrates both the single-atom model 
and the symmetric-sharing model, we need to find a good 
approximation for the 21 parameters in F. As insufficient 
theoretical information has been published on the frag- 
mentation of highly- charged N2 , we need to find the pa- 
rameters with the help of experimental data. The quest 
for the right parameters is further hindered by the fact 
that the nominal pulse energies specified differ from the 
actually available energies^ The discrepancy between 
the nominal and the actual value of the pulse energy at 
the sample in the experiment is caused by losses in the 
beam line transmission which are not particularly well 
determined. In Ref. 0, it is specified that only 15 %- 
35 % of the nominal pulse energy arrive at the sample. 
In order to determine the fragmentation constants, we fit 
the theoretical ion yields from the fragmentation-matrix 
model [Eq. ([5])] to the experimental ion yields of N2 by 
minimizing the modulus of the difference between the av- 
erage charge states © with respect to the fragmentation 
constants at the four measured x-ray pulse durations, 

1. e., for 4, 7, 80, and 280 fs^ Thereby, we adjust the 
pulse energies such that the constant for the fragmenta- 
tion of [Eq. ([TO])]. /{ ,2} = 0.74, i.e., the experimen- 
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tal value for fragmentation after Auger decay of a SCH is 
recovered. This procedure yields the fragmentation con- 
stants and the fractions of the nominal pulse energy that 
are listed in Table Q] Note that this approach assumes 
that the molecular charge state {0, 2} is predominantly 
produced by SCH decay, i.e., the competing processes of 
double valence ionization by x rays is neglected because 
of tiny cross sections. 

A number of insights about the experiment can be 
gained by inspecting TableHl Examining Ep, we find that 
the actually available pulse energy in the experiment is 
with one exception in the range 25-31% of the nomi- 
nal pulse energy Only the pulse energy for 7 fs pulses is 
around 16% of the nominal pulse energy of 0.26 mJ. This 
drop can be rationalized by taking 28% of 0.15 mJ — the 
nominal pulse energy specified for the case of 4 fs pulses — 
which leads to the same actual pulse energy as 16% 
of 0.26 mJ. The value 28% is consistent with the per- 
centages of the other pulse durations and indicates that 
potentially there was a drop in the actual output of 
LCLS when going to so short pulses. The range 25- 
31% is higher than the value of ~23% specified in 
Ref. H for 1050 eV photon energy. The peak x-ray in- 
tensities follow to 3.1 x 10 17 W/cm 2 , 1.7 x 10 17 W/cm 2 , 
2.6 x 10 16 W/cm 2 , and 9.0 x 10 15 W/cm 2 for the nominal 
pulse durations 4, 7, 80, and 280 fs, respectively. 

The change of molecular fragmentation with pulse pa- 
rameters can be read off Table HI Predominantly, /{0.2} is 
determined by the absorption of one x ray, /{o,3} an d 
/jo, 4} are from two-x-ray-absorption, /{2,4}i /{o,5}j an d 
/{o,6} are from three-x-ray-absorption, and /{2,5} is from 
four-x-ray-absorption. Apart from /{2,4} an d /{o,6} the 
fragmentation constants overall increase with decreasing 
pulse duration. This indicates that symmetric sharing 
becomes very important for the shortest pulses while 
molecular fragmentation may occur during the longer 
pulse prior absorption of further x rays. Due to very 
low probabilities, the values for /{2,4}, /{2.5}j an d /{0.6} 
for 4fs pulses are not important. As /{2,4} an d /{o,6} are 
basically unchanged for all pulse durations, the associ- 
ated fragmentation process is ascribed to the core-hole 
lifetimes which are independent of the pulse duration. 

The merits of the fragmentation matrix model are ex- 
hibited in Figs. H] and [5] where we plot q and the ion 
yields, respectively, from the model alongside the exper- 
imental data. The agreement of our theoretical q with 
the experimental q is very good, however, despite fit- 
ting the fragmentation matrix at all points, the two q do 
not coincide. This has several reasons. First, the fit is 
not arbitrary but physically motivated by the molecular 
breakup process relying on the probabilities to find the 
atom in various charge states. The atomic rate equations 
of Sec. [IlT] were designed carefully, yet only a one-electron 
approximation has been used and particularly shake off 
is not treated.— ^ Other works on neon atoms^ show a 
far less than perfect agreement of atomic ionyields from 
experiment with theory. In the light of Refs. and|^, our 
theoretical ion yields in Fig. [5] appear almost too close to 
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FIG. 4. (Color online) Average charge state q from N2 subject 
to LCLS x-ray pulses of varying FWHM duration. We show q 
deduced from the experimental ion yields of N2 [Fig. [5] for the 
nominal pulse durations 4, 7, 80, and 280 fs, alongside our cal- 
culations; the points are connected by interpolation. Specif- 
ically, we plot the experimental q as black circles, q from a 
single N atom as red squares [Sec. IIII| , q from the symmetric- 
sharing model as green diamonds [Sec. lIV A] , and q from the 
fragmentation-matrix model [Sec. lIV B] as blue triangles. See 
Table [J for further LCLS pulse parameters. 



the experimental ones. 

The finding that q from the fragmentation-matrix 
model is lower than the experimental q for longer pulses 
is ascribed to the fact that no fragmentation during the 
pulse is treated. Namely, the fragmentation time of N?, + 
is2Si ~400 fs and thus shorter than the longest, 280 fs, 
LCLS pulses. Hence molecular fragmentation of N? + 
may occur before the x-ray pulse is over. When N 2 + 
fragments into N + + N + , absorption of another x ray 
by the atomic ions N + leads to an enhancement of the 
N 3+ yield; fragmentation into N + N 2+ increases the N 2+ 
and N 4+ yields. The fra gmentation-matrix model does 
not distinguish this case from breakup of the respec- 
tive molecular ions which cannot be treated by static 
ratios. A time- and molecular-state-dcpcndcnt fragmen- 
tation matrix should be able to fully account for the ob- 
served molecular ion yields. 

Finally, we used the nominal pulse duration in all cal- 
culations so far as was done in Ref. 0. However, the 
pulse duration was found to be significantly shorter in ex- 
periments^^ Namely, a nominal pulse duration of 80 fs 
was estimated actually to be 20-40 fs in Ref. H and mea- 
sured to be 40 fs in Ref. |63|; the nominal pulse duration 
of 300 fs pulses was found to be 120 fsjSS Yet the experi- 
mental conditions in the two studies Refs. [H and |H differ 
somewhat from the one in Ref. 0- Fortunately, the depen- 
dence of our theoretical predictions on the pulse duration 
turns out to be quite low^ 
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FIG. 5. (Color) Ion yields Yj of N2 subject to LCLS x-ray pulses of varying duration: (a) 4fs, (b) 7fs, (c) 80 fs, and (d) 280 fs. 
Experimental Yj are given by the blue bars and theoretical Yj from the fragmentation-matrix model [Sec. IIV B| by the red 
bars. See Table [J for further LCLS pulse parameters. 



V. EXPLANATION OF THE OBSERVED 
PHENOMENA IN TERMS OF ELEMENTAL 
MOLECULAR PROCESSES 

In the previous Sees. IIII1 IIV Al and IIV B[ we inves- 
tigated phenomcnological models that are based on the 
solution of the atomic rate equations to describe the in- 
teraction of N 2 with x rays where we used the average 
charge state q [Eq. ©] as a figure of merit to assess 
the quality of our models. Here we discuss the problem 
from the perspective of the processes entering the atomic 
rate equations to unravel the elemental molecular pro- 
cesses that give rise to the behavior of the experimental q 
in Fig. |H In Refs. 1 and l3ll a molecular-rate-equation 
model has been discussed based on elemental molecular 
processes which differs in various aspects from the fol- 
lowing analysis. 

To describe the elemental molecular processes of N2 in 
intense x rays, we need to find out, first, which molecu- 
lar electronic configurations prevail in the course of the 
interaction with the x-ray pulse. As the atomic ioniza- 
tion cross sections for core electrons are orders of mag- 
nitude larger than cross sections for valence electrons, 
predominantly core electrons are ionized. The probabil- 
ities to find N2 in its ground state or with a SCH, a 



tsDCH, or a ssDCH can be deduced from the solution 
of the atomic rate equations for an N atom via P223(t) 2 , 

2Pi23(i)i 5 223(*), P123MP123W, and P023WP223W, re- 
spectively. Furthermore, a sequence of photoioniza- 
tion (P) with a following Auger decay (A) and another 
photoionization (PAP) is a competing process to DCHs 
for two-x-ray-photon absorption ! 11 ' 28 These probabilities 
are shown in Fig. [fj] for a short-pulse (4 fs) and a long- 
pulse (280 fs) scenario. It is revealed clearly that for 
the long pulse [Fig. [BJd], we have a contribution from 
almost exclusively SCHs and PAPs with an almost van- 
ishing probability to find DCHs. On the contrary, for 
the short pulse, PAPs are less important due to the fact 
that the Auger decay time of a SCH2£ of 6.7 fs is longer 
than the pulse duration. Therefore, DCHs are more rel- 
evant for the short pulse [Fig. [6^,] where tsDCHs make 
a larger contribution than ssDCHs. This is due to the 
bigger cross sections for the formation of tsDCHs com- 
pared with ssDCHs and the longer lifetime of the for- 
mer with respect to the latter. Namely, for an N atom 
in x-ray light at HOOcV photon energy, the cross sec- 
tion for SCH formation is 56kbarn and it is 34kbarn 
for DCH formation. The SCH formation cross section 
is, therefore, by a factor of 1.7 higher than the cross 
section for DCH formation in an N atom but ssDCHs 
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FIG. 6. (Color) The probabilities P(t) to find N2 during the x-ray pulse in its ground state is given by the black lines; the 
red lines are the probability for a SCH; the green and blue lines stand for the probability to find a tsDCH and a ssDCH, 
respectively; the PAPs are indicated by magenta lines. Probabilities are determined for (a) a Gaussian short pulse (J7J of a 
FWHM duration of 4fs and (b) a Gaussian long pulse of a FWHM duration of 280 fs. See Table [J for further LCLS pulse 
parameters. 



in N2, nevertheless, make a noticeable contribution for 
short pulses, i.e., higher x-ray intensities, as is revealed 
in Fig. ink- The lifetime of a tsDCH is half the lifetime 
of a SCH, i.e., 6.7fs/2 = 3.4 fs, versus 2.2 fs for a ssDCH, 
respectively [Tables IV and IX of Ref. [H]. The conclu- 
sions drawn from Fig. |S] are corroborated by the ratios 
in Table |TT] to find SCH decay, first decay of DCHs, and 
PAPAs after volume integration ©. We also show the 
impact of the SASE pulse structured— on the ratios 
using the partial coherence method to generate SASE 
pulse o 25 ' 64 ' 65 with 8 cV bandwidth and a Gaussian tem- 
poral envelope that has the same FWHM after squaring 
as the Gaussian pulse used to compute the reference ra- 
tios and make Fig. [51 Inspecting Table HH we find that 
the impact of the SASE pulse structure on the resulting 
ratios is small as has also been found in previous work»ii 
Another reason why SCH, ssDCH, tsDCH, and 
PAP channels have a particular importance is the fact 
that Nj" and N 2 , -1- are metastable with respect to dissocia- 
tion^^ i.e., the molecular ion remains essentially intact 
until it has charged up to N 3 , 4 " and higher.— Therefore, 
these channels dominate the first stages of the interac- 
tion of N2 with x rays and thus have a crucial impact on 
fragmentation products and such on molecular ion yields. 
Similarly, for very high intensities, also N2 with triple 
and quadruple core holes may occur. However, for our 
parameters, Fig. [5] reveals that we need to focus only on 
a very limit number of initial molecular channels. Based 
on the discussion of Sees. IIIIl IIV A[ and IIV Bl we deduce 
that one needs to consider following elemental molecular 
processes: 

First, for long pulses with moderate x-ray photon flux, 
where the probability to absorb two x rays in a few- 
femtosecond time interval is low, N 2 , -1- frequently frag- 
ments into two atomic ions prior to absorbing further 
x rays. However, for a higher x-ray flux this probability 
is not low and differences may arise due to further ab- 



sorption of x rays prior to molecular dissociation. This 
introduces a time scale of ^100 fs that determines the 
rate of fragmentation of the dication N 2+ [Ref. H3| • 

Second, for high enough x-ray intensities, further pho- 
toionization and Auger decay of the molecular dica- 
tion N 2+ prior fragmentation forms molecular trica- 
tions N| + . No bound states were found for N?, + ions 
in the midst of many Coulomb repulsive potential sur- 
facesi 59 ' 66 Therefore, one may assume that those N 3 , -1- 
without core holes fragment into N 2+ and N + which is 
the channel that is indicated most prominently by the 
potential energy surfaces. There is only a small amount 
of N 3 , -1- produced without core holes in the present case 
due to small valence ionization cross sections for x rays. 
Although N 3 ," 1- is not metastable, dissociation is not im- 
mediate and further x rays may be absorbed prior dis- 
sociation. Specifically, N?, 4 " with core holes should still 
be considered as a molecular entity which is only broken 
up after Auger decay because of the faster time scale of 
Auger decay compared with molecular fragmentation. 

Third, the sequence of two absorptions of an x rays 
with subsequent Auger decays (PAPAs) or, likewise, the 
two Auger decays of a DCH leads to ■ There is hardly 
any difference between the two situations concerning the 
way molecular fragmentation takes place. A SCH decay 
leads to N 2 ," 1 " which is metastable with a lifetime that is 
large compared with inner-shell hole lifetimes! 23 ' 28 ' 30 The 
absorption of a second x ray then leads to N| + which is 
not metastabl e 59 ' 66 and the molecular ion starts to disso- 
ciate, however, on a much slower time scale than Auger 
decay. The Auger decay time scale is independent of 
the x-ray pulse duration. As shown in Table [Hi PA- 
PAs are comparable with the first decay of ssDCHs in 
the 4fs case and PAPAs thus are an important channel 
for two-x-ray-photon absorption. Furthermore, only a 
DCH is metastable with respect to dissociation; as soon 
as the first Auger decay of a DCH occurs, N?, + is formed 
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rx [fs] 


Pulse 


SCH decay 


ssDCH decay 


tsDCH decay 


PAPA 


4 


1 SASE 


0.87 


0.049 


0.079 


0.052 




10 SASE 


0.88 


0.046 


0.073 


0.045 




100 SASE 


0.87 


0.050 


0.080 


0.045 




1000 SASE 


0.87 


0.050 


0.081 


0.046 




Gaussian 


0.89 


0.041 


0.065 


0.035 


280 


1 SASE 


0.99 


0.0035 


0.0058 


0.28 




10 SASE 


0.99 


0.0045 


0.0074 


0.33 




100 SASE 


0.99 


0.0044 


0.0072 


0.33 




Gaussian 


0.99 


0.0043 


0.0072 


0.33 



TABLE II. Ratios to find SCH decay, first tsDCHs decay, the first ssDCHs decay, and a PAPA process after averaging over a 
varying numbers of SASE pulses and from a Gaussian pulse with a FWHM duration of rx. The ratios are probabilities from 
volume averaging Q which are renormalized to the sum of the probabilities for (the first) decay of SCHs, tsDCHs, and ssDCHs. 
See Table U for further LCLS pulse parameters. 



which starts to dissociate in straight analogy to PAP af- 
ter the absorption of the second x ray. The second Auger 
decay of a DCH is on the same time scale as for the 
PAPA process. As both processes describe very similar 
situations with respect to molecular fragmentation, they 
have the same fragmentation pattern. From our find- 
ings with the symmetric-sharing model [Sec. IIV A| and 
results from laser-dissociation experiments j 60 i 61 we may 
assume equal sharing of the four charges between the two 
nitrogen nuclei for the fragmentation of ions. This 
is corroborated by inspecting Fig. [5^i where no signifi- 
cant ion yield for N 3+ is found which would arise from 
a nonsymmctric sharing of charges for the fragmentation 
of Nj" 1 " ions resulting from DCHs or PAPAs. In our case, 
DCHs and PAPs affect results for the short and long 
pulses where the importance of DCHs is larger for the 
high x-ray flux of short pulses because the molecules are 
ionized faster. Fragmentation predominantly occurs after 
the short pulses are over. For long pulses, the x-ray flux 
is reduced leading to a slower rate of photoabsorption 
and frequently fragmentation of N| + occurs before fur- 
ther x rays are absorbed. 

Fourth, for even higher charge states, dissociation 
becomes more and more rapid but the involved time 
scale needs, nonetheless, be considered. Namely, for 
x-ray pulses of only a few femtoseconds, dissociation does 
not play a substantial role and thus the assumptions un- 
derlying the fragmentation-matrix model [Sec. IIV Bj arc 
well fulfilled. In this case, also for highly-charged N2, 
symmetric sharing should be the predominant mode 
of fragmentation^ For long x-ray pulses, breakup fre- 
quently occurs before high molecular charge states are 
reached and subsequent x-ray absorption occurs by inde- 
pendent ions. 

The above analysis suggests that molecular effects and 
molecular fragmentation of N2 in intense x rays are cen- 
tered around these four steps. The described elemen- 
tal molecular processes should characterize the most im- 
portant aspects of the phenomenology of the interaction 



as all insights gained from the analysis of the models of 
Sees. IIIII IIV Al and IIV Bl are incorporated. 



VI. CONCLUSION 

In this paper, we discuss theoretically the interaction 
of intense, short-pulse x rays from the novel x-ray free 
electron laser Linac Coherent Light Source (LCLS) with 
N2 molecules. We set out from the computation of the 
energies, decay widths, and one-x-ray-photon absorption 
cross sections of a single N atom. These data are input 
to describe the interaction an N atom with x-ray pulses 
in terms of a rate-equation model. Ion yields and the 
average charge state are computed and compared with 
experimental data for N2. The agreement of theoretical 
data for an N atom and experimental data for N2 is not 
satisfactory which we ascribe to a fragmentation of N2 
after ionization causing a redistribution of charges over 
both atoms in the molecule. To describe the molecu- 
lar degrees of freedom, we devise two phcnomcnological 
models. First, we compute a lower limit to the average 
charge state assuming always equal sharing of charges be- 
tween both atoms in N2. Second, we use a fragmentation 
matrix to treat molecular dissociation with static ratios 
to relate atomic and molecular ion yields. The mod- 
els clearly indicates the relevance of the redistribution 
of charge and a different weighting of processes on short 
and long time scales for similar nominal x-ray pulse en- 
ergy. The formation of high charge states at short pulse 
durations is suppressed by a redistribution of charge in 
the multiply-ionized molecule giving rise to the observed 
frustrated absorption^ 

Pertinent to our atomic rate equations of Sec. IIIII how- 
ever, is — as was shown by Young et al£ for neon atoms — 
that the deviation of the theoretical ion yields from the 
experimental ones in the intense x rays of LCLS becomes 
more and more pronounced when the photon energy is of 
the order of 2000 eV. The reason of this deviation is the 
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fact that in Ref. H, the rate equations for a neon atom 
contain only all possible one-photon absorption and flu- 
orescence and Auger decay channels as is also the case 
for our rate-equation model of an N atom [Sec. IIII| . How- 
ever, at such high photon energies, other channels involv- 
ing many-electron effects like shake off (for neon atoms 
in Ref. M) may become important. As in the present 
study of N2, the photon energy of 1100 cV also is high 
above the SCH-ionization threshold of an N atom, a sim- 
ilar deviation, as was found for neon,— ^ may be relevant. 
However, we do not investigate this issue due to a lack of 
experimental data on single N atoms in intense x rays. 

Although our theoretical analysis is centered on N2 , we 
believe that many of our conclusions are generally true. 
The time scale of the involved processes is dictated by the 
formation and decay of SCHs and DCHs. Higher num- 
bers of core holes than one DCH per N 2 are of little im- 
portance at the present x-ray intensities. Molecular frag- 
mentation depends largely on the pacing of x-ray absorp- 
tion and, related, the abundance of SCHs versus DCHs. 
For x-ray irradiation durations of only a few femtosec- 
onds, nuclear dynamics — specifically due to Coulomb ex- 
plosion of the charged fragments — does not play a role 
during the x-ray pulse and we may assume fixed nuclei 
in the equilibrium geometry of the molecule. However, 
the nuclear wavepacket generated leads to a subsequent 
fragmentation of the molecule. Instead, for longer pulse 
durations, nuclear expansion during the x-ray pulse be- 
comes a crucial issue. A molecule with a SCH may be 
metastable with respect to dissociation and may even 
remain metastable after Auger decay has taken place. 
Clearly, this is crucial for the understanding of the ob- 
served ion yields. Hence molecular effects and molecu- 
lar fragmentation are centered around these initial steps. 
For a general description of molecules in intense x rays, 
one needs to account for the complicated joint problem of 
x-ray absorption, decay processes, and nuclear dynamics. 
We make first attempts towards a full treatment of the 
problem. 

Based on our studies, many fascinating possibilities 
for future research open up. An investigation with ab 
initio methods of molecular physics^ of the absorption 
of x rays, the decay of inner-shell holes and the result- 
ing nuclear dynamics poses a challenge. Nonetheless, we 
consider it of high importance to find an ab initio descrip- 
tion to gain deeper insights and to establish a theoretical 
foundation for molecules in intense x rays. The formation 
of DCHs gives one detailed information about the struc- 
ture of a molecule; it was found to be more sensitive than 
SCH spectroscop y 12 i 13 ' 15 i 16 The time-delay between the 
two core-hole ionization events in DCH production needs 
to be short enough, such that the first vacancy has not 
decayed. Then the second photoelectron carries detailed 
information about the dicationic states of the molecule. 
This gives rise to x-ray two-photon photoelectron spec- 
troscopy (XTTPS)f^ Furthermore, we have a valence 
electron density in molecules that may differ appreciably 
from the density produced by the noninteracting atoms 



in the molecular geometry. This may in some cases have a 
noticeable influence on the rate of intraatomic electronic 
decay (IAED)jSI if valence electrons are involved. Ad- 
ditionally, the close packing of the atoms in a molecule 
gives rise to interatomic electronic decay processes for 
inner- valence holes, namely, interatomic Coulombic de- 
cay (ICD) and electron transfer mediated decay (ETMD) 
(Ref. [6?] and References therein) . The relevance of these 
additional channels needs to be assessed in future work. 
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